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1.1  Nonlinear  Light  Mixing  Spectroscopy.  Y.  Prior,  A.  Bogdan,  M. 

Dagenais,  and  N.  Bloembergen,  Contract  N00014-75-C-0648;  Research 
Unit  5. 

The  work  on  nonlinear  four-wave  mixing  has  been  continued.  The 
emphasis  has  been  shifted  to  four-wave  light  mixing  in  sodium  vapor. 

Several  groups  working  on  the  theory  of  the  nonlinear  susceptibility,  * 

predict  slightly  different  results  for  the  lineshapes  of  CARS  (Coherent 
Anti-Stokes  Raman  Scattering)  and  CSRS  (Coherent  Stokes  Raman  Scattering) 

signals.  In  the  previous  annual  report  the  general  theoretical  expression 

1 2 
obtained  by  us  was  mentioned.  This  theoretical  work  was  also  reported 

at  the  Sixth  International  Conference  on  Raman  Spectroscopy,  Bangalore, 

India,  September  1978. 
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An  experimental  test  of  some  of  the  damped  resonance  terms,  about 

which  some  controversy  still  exists,  has  been  undertaken.  Atomic  sodium 

was  chosen  because  the  energy  levels  are  sharp  and  well  defined,  and  they 

are  connected  by  known  matrix  elements.  Furthermore,  the  damping  can  be 

varied  in  a controlled  manner  by  the  partial  pressure  of  a noble  buffer 

gas.  The  experiments  avoid  any  one-photon  resonance  from  the  populated 

3S  ground  state  of  Na,  but  the  frequencies  of  two  incident  dye  laser  beams 

are  at  or  near  resonance  between  the  excited  levels  of  the  3P..  . . doublet 

1/  2 1 3/  2 

and  the  4D  level.  The  beams  are  generated  in  a double  frequency  dye  laser 
of  a novel  design,  which  has  been  published.3  These  experiments  are 
difficult  to  perform  due  to  competing  effects  of  sodium  molecules  Na2» 
Four-wave  mixing  signals  have  been  detected  and  analyzed,  but  no  conclusive 
results  have  been  reached  to  date.  It  appears  that  a very  small  (one  part 
in  a million),  but  nevertheless  real,  population  exists  in  the  3P  atomic 
state.  Presumably  this  level  is  reached  following  molecular  excitation  of 
Na2«  Further  work  to  unravel  the  true  four -wave  mixing  by  the  ground 

state  population  is  in  progress. 

References 

1.  N.  Bloembergen , H.  Lotem,  and  R.T.  Lynch,  Jr.,  "Lineshapes  in  Coherent 
Resonant  Raman  Scattering,"  Indian  Journal  of  Pure  and  Applied  Physics, 
16,  151  (1978) . 

2.  N.  Bloembergen,  "Stimulated  Raman  Effect  and  Coherent  Raman  Scattering," 
Proc.  of  the  Sixth  International  Conference  on  Raman  Spectroscopy, 
edited  by  E.D.  Schmid,  R.S.  Krishnan,  W.  Kiefer,  and  H.W.  Schrotter 
(Heyden,  Philadelphia,  1978),  p.  335. 

3.  Y.  Prior,  "Double  Frequency  Narrow  Band  Grazing  Incidence  Pulsed  Dye 
Laser,"  Eev.  Sci.  Inst.,  50,  259  (1979). 
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1.2  Two-Photon  Absorption  in  Solids  and  Liquids.  Y.  Prior  and  H.  Vogt, 
Contract  N00014-75-C-0648;  Research  Unit  5. 

The  work  initiated  -y  Lotem  and  de  Araujo1  to  develop  accurate  methods 
for  measuring  two-photon  absorption  cross  sections  has  been  continued.  The 
method  of  measuring  TPA  absolute  cross  sections  by  means  of  calibration 
relative  to  known  Raman  cross  sections  was  developed  further  to  enable 
internal  calibration  in  Raman  active  liquids.  The  loss  coefficient  at  a>2 
is  proportional  to  the  pump  intensity  I (u^)  due  to  the  two-photon 
absorption  process  at  + U)^  ln  a crYstal  which  is  immersed  in  a fluid, 
e.g.,  benzene.  The  frequencies  ai  and  u are  so  chosen  that  w - m 

X A 1 4 

corresponds  to  a Raman  resonance  with  a known  Raman  cross  section  in 
benzene.  The  loss  due  to  two-photon  absorption  may  be  compared  to,  and 
compensated  by,  the  Raman  gain.  Both  processes  are  proportional  to  I(w^). 
The  method  is  accurate  because  it  is  independent  of  mode  structure,  self- 
focusing  and  other  distortions  of  the  light  beams.  The  measurements  were 
extended  to  the  uv  part  of  the  spectrum,  and  TPA  coefficients  were  obtained 
for  benzene  and  several  alkali  halide  crystals  at  excitation  energy  of 
6.7  eV. 

The  results  were  reported  at  the  1978  annual  meeting  of  the  Optical 

2 

Society  of  America  in  San  Francisco,  and  are  discussed  in  detail  in  a 

paper  accepted  for  publication  in  Phys.  Rev.  B.3  The  present  method 

utilizes  ns  pulses  derived  from  a uv  dye  laser  pumped  by  the  second 

harmonic  of  a ruby  laser.  Our  results  for  the  TPA  cross  sections  in  the 

uv  were  compared  to  those  obtained  in  this  laboratory  previously  with  ps 
4 

pulses.  There  is  good  agreement  between  the  two  independent  methods. 


W. 
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Phy8.  Rev.  B18,  30  (1978) . 
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Cross  Sections,"  J.  Opt.  Soc.  of  Amer.  , 68,  1423  (1978). 

3.  Y.  Prior  and  H.  Vogt,  "Measurement  of  uv  Two-Photon  Absorption  Relative 
to  Known  Raman  Cross  Sections,"  to  be  published  in  Phys.  Rev.  B,  June 
1979. 

4.  P.  Liu,  W.L.  Smith,  H.  Lotem,  J.H.  Bechtel,  N.  Bloembergen,  and  R.S. 
Adhav,  "Absolute  Two-Photon  Absorption  Coefficients  at  355  and  266  nm," 
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1.3  Two-Photon  Absorption  and  Optical  Damage  with  Picosecond  Light  Pulses. 
P.L.  Liu,  R.  Yen,  and  N.  Bloembergen,  Contract  N00014-75-C-0648  and 
O.N.R.  Grant  N60530-78-C-0041;  Research  Unit  6. 

Measurement  of  the  dielectric  breakdown  threshold  in  diamond  was 

carried  out  at  1.06  ym.  The  result  suggests  that  the  breakdown  mechanism 

in  diamond  is  similar  to  that  in  ionic  crystals.  The  upper  limit  of  the 

two-photon  absorption  coefficient  of  diamond  was  also  determined  at  532  nm. 

A cumulative  damage  effect  at  532  nm  was  observed.  It  was  a result  of  the 

accumulation  of  the  color  centers  generated  through  two-photon  absorption 

in  laser  pulses.  A paper  based  on  these  measurements  has  been  published.* 

Measurements  on  the  two-photon  absorption  coefficients  of  ultraviolet 

window  and  coating  materials  were  carried  out.  The  two-photon  absorption 

was  determined  in  fused  silica,  SrF2,  BaF2 , MgF2,  and  LiYF^ . The  data  are 

relevant  for  high  power  uv  lasers.  A paper  presenting  the  design  of  an 

automatic  data  acquisition  system  along  with  the  two-photon  absorption 

2 

coefficients  was  accepted  for  publication  in  Applied  Optics. 
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1.  P.  Liu,  R.  Yen,  and  N.  Bloembergen,  "Dielectric  Breakdown  Threshold, 
Two-Photon  Absorption,  and  Other  Optical  Damage  Mechanisms  in  Diamond," 
IEEE  J.  Quantum  Electron.  QE14,  574  (1978) . 

2.  P.  Liu,  R.  Yen,  and  N.  Bloembergen,  "Two-Photon  Absorption  Coefficients 
in  uv  Window  and  Coating  Materials,"  to  be  published  in  the  April  1,  1979 
issue  of  Applied  Optics. 

This  material  is  also  contained  in  the  Final  Technical  Report,  Feb. 

1979,  of  ONR  Grant  N60530-78-C-0041. 


1.4  Ultra  Fast  Quenching  and  Photoelectric  Emission  Induced  by  Picosecond 


Light  Pulses.  P.L.  Liu,  R.  Yen,  and  N.  Bloembergen,  Contract 
N00014-7 5-C-0648  and  ONR  Grant  N60530-78-C-0041;  Research  Unit  6. 


Picosecond  laser  pulses  are  used  to  induce  heating  and  cooling  with 
13 

a rate  of  the  order  of  10  °C/sec  on  the  surface  of  a single  crystal  of 
silicon.  Preliminary  results  demonstrate  a crystalline  to  amorphous  trans- 
formation induced  by  the  short  pulse.  This  ultra  fast  quenching  process 
may  lead  to  new  phenomena  on  the  kinetics  of  phase  transitions.  It  may 
also  find  application  in  optical  memory  switching.  Detailed  studies  of 
the  morphology  of  the  phase  transformation  are  planned  for  the  next 
reporting  period. 

The  frequency  and  polarization  dependence  of  the  multiphoton  photo- 
electric emission  are  under  study.  A tunable  parametric  oscillator  has 
been  designed  and  is  under  construction.  In  order  to  pump  the  parametric 
oscillator,  we  developed  a mode-locked  unstable  Nd : YAG  oscillator.  The 
output  pulse  has  a duration  of  28  picosecond  and  an  energy  content  of  8 mJ 
in  each  train.  A paper  presenting  the  result  of  the  unstable  oscillator 
was  accepted  for  publication.1 
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Reference 

1.  P.  Liu  and  R.  Yen,  "Passively  Mode-Locked  Nd : YAG  Unstable  Oscillator," 
to  be  published  in  the  March  1,  1979  issue  of  Applied  Optics. 


1.5  Intramolecular  Vibrational  Relaxation  in  Polyatomic  Molecules. 


I.  Burak,  H.S.  Kwok,  R.  Sharp,  J.  Tsao,  and  E.  Yablonovitch , Contracts 
N00014-7 5-C-0648  and  N00014-78-C-0531 ; Research  Unit  7. 


Probably  one  of  the  most  important  problems  in  chemical  physics  is 
the  concept  of  intramolecular  vibrational  relaxation.  How  fast  will  energy 
redistribute  from  one  vibrational  mode  to  the  other  vibrational  modes  of 
the  same  molecule?  This  question  has  been  fraught  with  controversy  for 
decades,  but  we  now  feel  that  the  modern  developments  in  laser  technology 
will  begin  to  answer  those  questions  in  the  next  few  years. 

The  problem  is  important,  for  quantum  electronics,  statistical 
mechanics  and  physical  chemistry.  It  is  also  central  to  infrared  multi- 
photon excitation.  In  the  past  year  we  completed  our  initial  experiments 
in  this  field. ^ Using  a time  resolved  infrared-infrared  double  resonance 
technique,  we  have  placed  an  upper  limit  of  30  psec  for  the  energy  re- 
distribution time  as  the  energy  relaxes  out  of  the  initially  excited 
mode  of  SF^ . A standard  picosecond  pump  and  probe  technique  was  used. 

Dr.  H.S.  Kwok  received  the  Ph.D.  for  this  experiment  and  has  taken  a post- 
doctoral position  at  Lawrence  Berkeley  Laboratory.  Mr.  Sharp  is  continuing 
those  experiments,  with  the  initial  objective  of  taking  a full  double 
resonance  spectrum  and  improving  the  signal-to-noise  ratio  by  computerizing 
the  data  collection. 
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Dr.  I.  Burak  and  Jeff  Tsao  are  taking  a different  approach,  and  have 
set  up  an  infrared-visible  double  resonance  system  based  on  our  picosecond 
CX>2  laser  pulses  and  a ruby  pumped  dye  laser  system.  They  have  synchronized 
the  two  laser  systems  within  25  nsec  at  present  and  are  already  conducting 
experiments  on  NC>2 . 

In  the  course  of  this  work  we  have  conceived  an  ingenious  digital 

feedback  technique  for  the  single  axial  mode  stabilization  of  CC>2  TEA 

2 

lasers.  This  work  will  be  published  shortly  and  is  expected  to  become 
the  standard  axial  mode  stabilization  technique  in  our  laboratory. 
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1.6  Multiphoton  Vibrational  Excitation.  J.  G.  Black  and  E.  Yablonovitch, 
Contracts  N00014-75-C-0648  and  N00014-78-C-0531;  Research  Unit  7. 

One  of  the  major  questions  in  infrared  laser  induced  chemical 
reactions  is : What  is  the  vibrational  energy  distribution  produced  by 
multiphoton  excitation?  A simplified  theory  and  experimental  results1 
t for  SF&  have  been  published  with  the  following  conclusion:  The  energy 

distribution  is  bimodal  with  some  fraction  of  the  molecules  bottlenecked 
among  the  discrete  levels  and  the  balance  excited  in  roughly  a thermal 
distribution. 

This  model  for  multiphoton  excitation,  in  which  energy  fluence 
rather  than  peak  power  plays  the  key  role,  was  described  in  a recent  review 
article. 2 
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During  these  experiments  we  discovered  a decomposition  formula  for 
Gaussian  spatially  averaged  experimental  data.3  The  true  nonlinear 
functional  dependence  of  a physical  parameter  on  local  light  intensity  can 
be  obtained.  This  mathematical  formula  is  invaluable  for  our  energy 
deposition  and  yield  experiments.  These  measurements  are  being  extended 
to  a series  of  molecules  of  different  sizes.  To  this  end,  the  data 
collection  system  has  been  computerized  and  tested.  Measurements  should 
be  completed  during  the  next  reporting  period. 
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Bloembergen,  "Collisionless  Multiphoton  Energy  Deposition  and 
Dissociation  in  SFg."  Phys.  Rev.  A,  February  (1979). 

3.  N.  Bloembergen  and  E.  Yablonovitch,  "Infrared  Laser  Induced  Unimolecular 
Reactions,"  Physios  Today,  31,  no.  5,  p.  23,  May  (1978). 

4.  P.  Kolodner,  H.S.  Kwok,  J.G.  Black,  and  E.  Yablonovitch,  "Exact 
Decomposition  of  a Gaussian -Averaged  Nonlinear  Function,"  Opt.  Lett., 

4,  38  (1979). 

1.7  Laser  Plasma  Interaction.  P.  Kolodner  and  E.  Yablonovitch,  Contracts 
N00014-75-C-0648  and  DOEU78-S-02-4631-A000;  Research  Unit  7. 

We  have  continued  to  study  the  resonant  acceleration  of  electrons 
in  the  interaction  of  a powerful  electromagnetic  wave  with  a plasma  density 
gradient.  This  effect,  which  was  first  experimentally  discovered  in  our 
laboratory,  plays  a major  role  in  the  absorption  of  laser  energy  by 
inertial  confinement  fusion  targets. 
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In  the  resonant  ids  sorption  (RA)  theory  of  laser-plasma  interaction, 

the  most  important  plasma  parameter  is  the  electron  density  scale  length 

L = n IVn  | 1.  Our  efforts  in  the  past  year  have  been  directed  towards 
6 © 

plasma  production  with  reproducible  L by  optical  breakdown  on  the  front 
of  a moving  shock  wave.  Such  a target  has  L~1  to  10  ym,  exactly  the 
range  of  interest  in  present  laser  interaction  research. 

To  this  end,  we  have  constructed  an  electrothermal  shock  tube  which 
produces  waves  of  Mach  number  ~5  in  H2  gas  at  ~50  torr  pressure.  A high- 
speed laser  schlieren-type  shock-front  detection  system  has  been  developed 
which  allows  the  motion  of  the  shockfront  to  be  "frozen"  with  spacial 
resolution  <100  ym.  Experiments  have  been  performed  in  which  high-power 
C02  laser  pulses,  triggered  by  this  detection  system,  irradiate  the 
"frozen"  shock  front.  The  spacial  density  profile  of  the  shock  front  has 
at  this  stage  manifested  itself  only  to  the  lowest  order;  i.e.,  when  the 
shock  front  is  far  from  the  focus  of  the  laser  pulse,  resonantly  accelerated 
electron  emission  is  greatly  reduced.  We  are  presenting  involved  in  attempts 
to  expose  more  directly  the  effects  of  the  neutral-gas  density  profile  in 
the  shock  front  on  the  laser-plasma  interaction. 
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II. 1 Defects  in  Liquid  Crystals.  S.A.  Asher,  J.  Collett,  and  P.S.  Pershan, 
Contracts  N00014-75-C-0648,  NSF  DMR77-24295  and  NSF  DMR76-22452; 
Research  Unit  4. 

The  work  described  last  year  on  defects  in  lipid  water  systems  has 

1 2 

been  continued  and  the  first  reports  on  this  topic  are  published.  ' 

Specific  features  of  the  proposed  parabolic  focal  conic  (PFC)  defect  have 
been  confirmed  experimentally  and  several  other  defects  have  been 
identified.  We  are  continuing  our  investigation  of  the  rheology  of  the 
PFC  defects  by  examining  their  formation  and  their  light  scattering 
properties  upon  mechanical  dilation  or  compression  of  monodomain  samples. 
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1.  S.A.  Asher  and  P.S.  Pershan,  "Parabolic  Focal  Conics  and  Polygonal 
Textures  in  Lipid  Liquid  Crystals,"  J.  de  Physique,  40,  161  (1979). 

2.  S.A.  Asher  and  P.S.  Pershan,  "Defect  Structures  in  Lyotropic  Lipid 
Liquid  Crystals,"  (to  be  presented  at  the  International  Colloid  and 
Surface  Science  Symposium,  Lolla,  MO,  June  1979)  . 


II. 2 Elastic  Properties  of  Liquid  Crystals.  M.  Fisch  ->nd  P.S.  Pershan, 

Contracts  N00014-75-C-0648,  NSF  DMR77-24295,  and  NSF  DMR76-22452; 

Research  Unit  4. 

Our  work  on  the  elastic  properties  of  the  smectic  phases  was  pre- 
viously based  on  the  Brillouin  technique.  Although  that  will  remain  an 
essential  part  of  our  effort,  we  have  investigated  the  feasibility  of 
studying  these  same  modes  at  lower  frequencies  using  a surface  wave 
technique.  We  have  learned  to  make  ~100  ym  thick  smectic  samples  with 
one  free  surface.  Alignment  and  optical  quality  of  these  samples  are 
better  than  any  others  that  have  previously  been  studied.  Theoretical 
considerations  imply  that  MHS  shear  acoustic  waves  in  the  bulk  of  these 
materials  form  underdamped  standing  waves  that  distort  the  free  surface 
and  thus  scatter  light  rather  strongly.  Preliminary  experiments 
demonstrate  that  the  theory  is  correct  and  we  are  currently  continuing 
to  develop  them  in  order  to  make  quantitative  measurements  of  the  low 
frequency  smectic  elastic  constants.  The  acoustic-optic  effects  in  this 
type  of  geometry  have  never  been  measured  and  we  will  evaluate  them  for 
applications  purposes  as  well  as  a tool  for  understanding  the  materials. 
Further,  these  results  together  with  Brillouin  studies  on  the  same 
material  will  give-  a clear  measure  of  the  importance  of  dispersive 
phenomenon  on  acoustic-optic  effects  in  smectic  liquid  crystals. 
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II. 3 Raman  Scattering  in  Liquid  Crystals.  S.A.  Asher  and  P.S.  Pershan, 

Contracts  N00014-75-C-0648,  NSF  DMR77-24295,  and  NSF  DMR76-22452; 

Research  Unit  4. 

The  construction  and  testing  of  the  PDP  11/03  computer  controlled 
Raman  spectrometer  has  been  completed.  Preliminary  measurements  on  lipid 
liquid  crystals  prepared  from  dipalmitoyl  phosphatidyl  choline  and  water 
in  the  Ux  phase  indicate  an  almost  isotropic  hydrocarbon  ordering.  How- 
ever, below  the  Lot  phase  the  molecular  ordering  appears  to  become  aniso- 
tropic. We  are  attempting  to  prepare  monodomain  samples  below  the  Lot 
phase  of  sufficient  optical  quality  for  polarized  Raman  measurements. 


II. 4 High  Frequency  Josephson  Phenomena.  D.A.  Weitz,  M.  Octavio,  W.J. 

Skocpol,  and  M.  Tinkham,  Contracts  N00014-75-C-0648,  N00014-76-C-0032, 
and  NSF  Grant  DMR76-11323;  Research  Unit  2. 

We  have  completed  systematic  measurements^"  of  the  ac  Josephson 

2 

effect  in  Nb  point  contacts  at  far-infrared  laser  wavelengths  between 

1.22  mm  and  70  pm.  Detailed  analysis  of  the  observed  current  widths  of 

the  constant-voltage  Josephson  steps  has  allowed  us  to  study  the  frequency 

dependence  of  the  strength  of  the  Josephson  effect  from  far  below  to  far 

3 4 

above  the  superconducting  energy  gap  frequency.  ' At  low  frequencies 
the  observed  step  widths  are  consistent  with  the  predictions  of  the  usual 
RSJ  model,  with  the  strength  of  the  Josephson  current  set  by  the  observed 
dc  critical  current.  Near  to  the  gap  frequency,  the  step  sizes  are  larger 
than  predicted  by  the  RSJ  model  in  a manner  qualitatively  consistent  with 
the  Riedel  peak  in  the  frequency  response.  Beyond  the  gap  frequency,  the 
size  of  the  steps  rolls  off  steadily  but  remains  easily  observable  at 


2.52  THz  (119  pm;  about  twice  the  gap  frequency).  At  the  highest  frequency 

studied  (4.3  THz  or  70  pm)  no  steps  were  observed,  but  the  primary 

difficulty  was  insufficient  laser  output  and  poor  coupling.  The  frequency 

dependence  generally  agrees  with  Werthamer's  tunneling  theory,6  but 

systematic  discrepancies  suggest  an  additional  rolloff  associated  with 

the  relaxation  time  of  the  order  parameter,  qualitatively  consistent  with 

predictions  of  the  time-dependent  Ginzburg-Landau  theory  applied  to  a 

metallic  constriction.6  The  data  definitely  excludes  the  need  for  the 

anomalously  long  relaxation  times  conjectured7  in  order  to  explain  data 
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from  indium  microbridges,  which  we  believe  were  affected  by  heating.  ' 
Certain  features  of  the  dc  I-V  curves  of  our  point  contacts, 
particularly  the  sharpness  of  the  structure  at  the  gap  voltage,  were 
found  to  correlate  well  with  good  high-frequency  performance  and  can  be 
used  as  reliable  predictors  of  high-frequency  behavior.  Other  features, 
such  as  the  "excess  supercurrent,"  were  shown  to  have  no  relation  to  the 
Josephson  effect.^"6 

We  have  also  investigated  the  noise  properties  of  our  junctions. 

The  high-frequency  noise  can  be  determined  from  the  rounding  of  the  steps 
induced  by  the  laser,  and  we  found  that  the  results  were  consistent  with 
Johnson  noise  enhanced  because  of  the  temperature  rise  of  the  point 

2 

contacts  under  irradiation.  Additional  noise  which  varies  like  1/f 
below  10  kHz  was  measured  directly,  and  proved  to  be  the  limiting  factor 


on  the  NEP  of  jucntions  used  as  incoherent  detectors  of  low-level  radiation 
when  biased  near  the  incipient  Josephson  step.  Typical  NEP  values  at 
604  GHz  (496  ym)  were  about  10  33  Watts/ (Hz) 1//2  with  a responsivity  of 
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II. 5 Laser-Induced  Disequilibrium  in  Superconducting  Films.  L.N.  Smith, 
A.D.  Smith,  W.J.  Skocpol,  and  M.  Tinkham,  Contract  N00014-75-C-0648, 
and  NSF  Grant  DMR76-11323;  Research  Unit  2. 

We  are  continuing  with  experiments  involving  optical  and  far- 
infrared  laser  illumination  of  superconducting  thin  films,  described  in 
greater  detail  last  year.  For  strong  optical  perturbation,  the  question 
being  investigated  is  whether  the  observed  smooth  transition  to  the 
normal  state  is  the  result  of  a spatial  instability  involving  breakup 
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into  domains  with  large  and  small  values  of  the  superconducting  energy 
gap. ^ Our  measurements  of  the  resistive  transition  have  proved  incon- 
clusive because  of  inhomogeneities  (phase-slip  centers)  are  also  caused 
by  the  strong  current  perturbation.  Most  of  our  recent  efforts  have  been 
devoted  to  developing  photolithographic  techniques  for  producing  small 
tunnel  junctions  which  are  directly  sensitive  to  the  energy  gap  and  which 
will  perturb  the  illuminated  film  less  strongly. 

We  are  also  continuing  our  efforts  to  measure  changes  in  the 
distribution  function  of  thin  films  illuminated  by  far  infrared  laser 
radiation  with  photon  energies  comparable  to  the  energy  gap.  As  a probe 
we  are  now  using  Al/I  Oxide/Pb-Bi  tunnel  junctions  which  are  sensitive  to 
the  distribution  function  over  a wider  range  of  energy  above  the  energy 
gap  of  aluminum.  We  have  carefully  characterized  the  equilibrium 
tunneling  curves  and  are  studying  the  distribution  function  with  weak 
optical  illumination.  We  have  not  yet  succeeded  in  absorbing  sufficient 
far-infrared  radiation  in  the  films  to  produce  interesting  nonthermal 
distributions  that  contain  structure  reflecting  the  photon  energy. 
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Anderson,  P.  Ketchian,  H.A.  Paesler,  and  W.  Paul.  Contracts 

N00014-7 5-C-0648 , NSF  DMR78-10014,  NSF  DMR76-15325,  NSF  DMR76-01111 

NSF  DMR77-24295  and  ARMY  DAAG29-77-G-0059;  Research  Unit  3. 

Amorphous  semiconducting  films  of  Ge,  Si,  Ge  H alloy.  Si  H 

a X X i 3C 

alloy,  GaAs,  and  (GaAs)^_x  alloy  have  been  produced  by  r.f.  sputtering 
of  polycrystalline  targets  by  argon  onto  substrates  of  Be,  7059  glass, 
crystalline  Si  and  Al  held  at  suitable  temperatures  between  o°C  and  400°C. 
Variation  of  the  partial  pressure  of  Ar  was  done  to  alter  the  conditions 
in  the  plasma  so  as  to  change  the  deposition  rate  and  the  surface  mobility 
of  the  depositing  Si  atoms.  Variation  of  the  substrate  temperature 
permitted  healing  of  the  disordered  network  without  inducing  crystalliza- 
tion. Variation  of  the  partial  pressure  of  H both  permitted  studies  of 
variable  composition  alloys  and  exploration  of  the  possible  different 
H-configurations  incorporated.  Addition  of  PH^,  B2H6'  N2  or  °2  the 
sputtering  gas  permitted  doping  the  semiconductor  n or  p-type,  the 
formation  of  n-p  junctions,  or  interfacing  high  and  low  resistivity 
material . 

The  films  are  characterized  chemically,  when  required,  by  electron 
microprobe  and  sputter-ion  analysis.  X-ray  examinations  are  conducted 
to  ensure  amorphicity,  the  density  is  determined  by  floating  the  sample 
in  a liquid  whose  density  may  be  varied  and  accurately  measured,  and  the 
thickness  is  found  using  a Sloan  Dektak  profiling  apparatus. 

The  films  produced  are  used  for  measurements  of  transport,  photo- 
transport, photoluminescence,  optical  absorption,  electron  spin  resonance, 
electroreflectance,  drift  mobility  and  field  effect.  A coherent 
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self-consistent  explanation  of  these  data  is  then  sought.  In  subsequent 
sections,  only  those  measurements  supported  predominantly  by  the  JSEP 
will  be  described. 


II. 7 Study  of  Doping,  Schottky  Barrier  Formation  and  Solar  Cell 

Efficiency  in  Amorphous  Si^  Alloy.  D.A.  Anderson,  G.  Moddel, 
and  W.  Paul,  Contracts  N00014-75-C-0648,  NSF  DMR76-01111,  NSF 
DMR77-24295,  and  DAAG29-77-G-0059;  Research  Unit  3. 


The  study  of  Schottky  barriers,  rectification  characteristics  and 
solar  cell  parameters  has  been  temporarily  postponed  in  favor  of  further 
determination  of  the  optimum  conditions  for  sample  preparation,  which  is 
described  in  Section  II. 8. 


II. 8 Study  of  Conditions  for  the  Optimum  Configuration  of  Pseudogat 


States  in  Amorphous  Silicon.  D.A.  Anderson,  R.  Collins,  E.C. 
Freeman,  G.  Moddel,  S.  Oguz,  M.A.  Paesler,  and  W.  Paul,  Contracts 
N00014-75-C-0648,  DAAG29-77-G-0059,  NSF  DMR76-15325,  and  NSF 
DMR78-10014 ; Research  Unit  3. 


During  this  report  period,  work  has  concentrated  on  further 
improvement  of  the  basic  material  so  as  to  optimize  the  functional 
dependence  (on  energy)  of  the  density  of  states  in  the  energy  gap.  In 
the  past,  optimization — whether  directed  at  producing  material  which  can 
be  doped  with  the  least  addition  of  dopant  atoms,  or  at  producing  material 
with  the  highest  photoconductivity  or  best  photoluminescence — has  been 
equated  to  the  maximum  reduction  of  defect-associated  states  in  the 


semiconductor  pseudogap.  It  now  appears  that  this  is  too  simplistic  a 
view,  as  will  be  described.  For  this  investigation,  the  substrate 
temperature  Tg,  the  H partial  pressure  pH»  the  argon  pressure  p^  and 
the  r.f.  power  have  all  been  varied. 

Evaluation  of  the  results  of  the  variation  of  these  deposition 
parameters  is  done  through  study  of  the  variations  in  conductivity,  thermo- 
power,  absorption  edge,  vibrational  absorption,  H-evolution,  photo- 
luminescence and  field-effect.  A detailed  discussion  of  these  results 
is  not  appropriate  here,  but  the  conclusions  can  be  given  in  summary 
form. 

Increasing  the  substrate  Tg  might  be  expected  to  heal  the  net- 
work, and  it  does,  but  high  T (>  300°C)  leads  also  to  less  incorporation 

s 

of  H for  a fixed  p . As  a result,  dangling  bond  defects  are  not 
adequately  compensated.  The  best  Tg  now  appears  to  be  the  highest 
temperature  which,  while  not  inducing  crystallization,  still  alloys  in- 
corporation of  H to  eliminate  dangling  bonds.  This  means  that  Tg  and 

p must  be  simultaneously  increased. 

H 

The  H-content  C is  not  a linear  function  of  the  H-partial 
H 

pressure  p . In  our  apparatus,  C„  increases  with  p for 
H H H 

p < ^.0  3 TOrr  and  C either  remains  constant  or  even  decreases  slightly 
H " H 

-3 

for  p >~10  Torr.  This  is  attributed  (as  was  reported  last  year)  to 
H " 

changes  in  the  nature  of  the  atomic  and  molecular  species  in  the  plasma 

when  the  gas  content  is  changed  with  r.f.  current  held  constant.  Films 

of  the  same  C produced  at  different  p have  been  found  to  possess 
H H 

quite  different  properties.  This  has  been  one  of  the  most  significant 
discoveries  of  this  report  period.  Films  prepared  at  low  p have  a 
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smaller  density  of  gap  states,  as  determined  by  measurements  of  hopping 
transport  through  them,  and  field  effect.  They  are  therefore  most 
suitable  for  shift  of  the  Fermi  level  by  P,  As  or  B doping,  and  they 
show  the  best  photoluminescence.  Films  of  the  same  C but  prepared  at 

high  p„  have  a much  higher  (order  of  magnitude)  density  of  gap  states, 

n 

as  revealed  by  hopping  transport  and  field  effect.  However,  interestingly, 
they  uniformly  give  a higher  photoresponse.  We  hypothesize  that  the 
distribution  of  states  in  the  gap,  determining  the  Fermi  level,  and  also 
trapping  and  release,  is  important  for  the  photoconductive  response. 
Obviously  this  has  direct  relevance  to  optimization  for  solar  cell 
applications. 

The  effect  of  varying  the  argon  pressure  has  also  been  studied, 
predicated  on  the  supposition  that  damage  to  the  growing  film  can  result 
from  very  energetic  Si  or  Ar  atoms/ions.  Increased  p^  reduces  the  mean 
free  paths  of  the  depositing  species  below  the  inter-electrode  distance 
and  so  reduces  their  kinetic  energy.  On  the  other  hand,  too  high  a p^ 
may  reduce  the  ions'  energy  to  thermal  values,  which  may  be  counter- 
productive, in  that  the  low  surface  mobility  can  lead  to  "fluffy"  porous 
films.  It  is  indeed  found  that  an  optimum  exists,  all  other 

conditions  being  held  constant.  Thus,  for  example,  the  photoresponse 
rises  then  saturates  around  pftr = 30  m Torr,  the  transport  remains  activated 
with  reduced  activation  energy  up  to  pftr  ~ 30  m Torr  but  above  this  shows 
an  increasing  curvature  with  temperature,  and  the  photoluminescence  also 
goes  through  a weak  maximum  in  this  pressure  range.  Furthermore  films 
made  at  p^r  > 30  m Torr  show  a tendency  to  oxidize  after  deposition 
suggesting  a large  surface  area  as  might  be  expected  for  an  extremely  porous 


II. 9 Use  of  the  Field  Effect  to  Determine  the  Pseudoqap  State  Density  in 

Amorphous  Si^^H^  Alloys.  D.A.  Anderson  and  W.  Paul,  Contracts 
DAAG-29— 77— G-0059  and  N00014— 75— C— 0648;  Research  Unit  3. 

During  this  report  period,  we  have  developed  an  apparatus  to  estimate 
the  density  of  states  in  the  gap  of  sputtered  a-Si^  alloys  through 
measurement  of  the  field  effect.  In  brief,  the  field  effect  technique  uses 
a configuration  similar  to  an  MOS  device  to  deposit  charge  in  the  surface 
layers  of  a semiconductor.  The  amount  of  charge  is  given  by  the  applied 
field  and  the  capacitance  of  an  insulating  dielectric.  The  charge  raises 
the  Fermi  level  at  the  surface  by  filling  the  states  at  the  Fermi  level. 

The  shift  in  the  Fermi  level,  in  turn,  is  found  from  the  change  in 
activation  energy  for  conduction  parallel  to  the  altered  surface  layer. 
Despite  assumptions  and  questionable  simplifications  in  the  analysis 
practiced  generally  in  the  field,  the  method  has  been  widely  used  to  give 
the  gap  density  of  states  in  a-Si^  alloys. 

Our  measurements  have  been  carried  out  on  alloys  of  different  H- 
content.  The  deduced  density  of  states  profile  depends  on  the  H-content, 
c„,  and  also  on  the  pressure  of  H gas,  p„,  used  to  obtain  a given  c„. 

The  samples  which  give  a high  gap  density  of  states  also  show  structure 
in  the  photoluminescence  spectrum  and  high  hopping  conductivity  which  is 
consistent  with  the  field  effect  data.  Thus  it  seems  clear  that  the 
addition  of  this  technique  will  help  our  overall  research. 


II. 10  Photoconductivity  in  Amorphous  Si^^H^  Alloys.  D.A.  Anderson, 

G.  Moddel,  and  W.  Paul,  Contracts  N00014-75-C-0648  and  DAAG-29-77- 
G-0059?  Research  Unit  3. 

We  have  extended  our  earlier  measurements  to  determine  the  spectral 

response,  the  temperature  dependence,  and  the  time  decay  characteristic  of 

photoconductivity  of  samples  of  a-s^i_xHx  all°Ys  prepared  under  different 

conditions  of  H partied,  pressure.  From  the  photoconductivity  magnitude 

and  the  decay  time  after  illumination  is  removed,  we  have  estimated  the 

-2  -3  2 

(trap-limited)  drift  mobility,  which  is  between  10  and  10  cm  /Volt-sec. 
for  electrons  at  room  temperature . There  are  significant  variations  of 
the  decay  time  with  illuminating  light  intensity,  whose  investigation 
continues . 

In  parallel  with  this  work,  a collaborative  effort  has  been  started 
with  Professor  John  Brebner  of  the  University  of  Montreal  to  determine  the 
drift  mobility  by  direct  time-of -flight  methods.  A report  of  the  pre- 
liminary success  of  this  work  will  be  given  by  members  of  the  two  groups 
at  the  next  meeting  of  the  Canadian  Physical  Society. 


Publications  of  the  Research  Reported  in  Sections  II. 6 to  11:10. 


1.  M.A.  Paesler  and  W.  Paul,  "Photoluminescence  of  r.f.  Sputtered 
Hydrogenated  Amorphous  Silicon,"  Bull.  Am.  Phys.  Soc. , 23,  247  (1978). 

2.  M.A.  Paesler  and  W.  Paul,  "Photoluminescence  in  Sputtered  Amorphous 
Silicon-Hydrogen  Alloys,"  submitted  to  Phil.  Mag. 

3.  S.  Oguz  and  M.A.  Paesler,  "Evolution  of  Hydrogen  from  Reactively 
Sputtered  Amorphous  Silicon,"  Bull.  Am.  Phys.  Soc.,  23,  247  (1978). 
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4.  W.  Paul  and  E.C.  Freeman,  "Infrared  Vibrational  Spectra  of  Si-H  Bonds 
in  a-Si,"  Bull.  Am.  Fhye.  Soo. , 23,  249  (1978). 

5.  E.C.  Freeman  and  H.  Paul,  "Infrared  Vibrational  Spectra  of  r.f. 
Sputtered  Hydrogenated  Amorphous  Silicon,"  Phy8.  Rev.  B18,  4288  (1978) 

6.  E.C.  Freeman  and  W.  Paul,  "Optical  Absorption  Edge  of  Sputtered  Hydro- 
genated Amorphous  Silicon,"  Bull.  Am.  Phye.  Soo.,  23_,  249  (1978). 

7.  M.A.  Paesler,  D.A . Anderson,  E.C.  Freeman,  G.  Moddel,  and  w.  Paul, 

"A  New  Development  in  the  Study  of  Amorphous  Si-H  Alloys:  The  Story 
of  0,"  Phye.  Rev.  Lett.,  41,  1492  (1978). 

8.  W.  Paul,  "Amorphous  Semiconductors,"  Proceedings  of  the  Third  Inter- 
national Conference  on  Physics  and  Contemporary  Needs,  Nathiagali, 
Pakistan,  to  be  published  by  Plenum  Press. 

9.  E.C.  Freeman  and  W.  Paul,  "Optical  Constants  of  r.f.  Sputtered  Hydro- 
genated Amorphous  Silicon,"  submitted  to  Phys.  Rev. 

10.  D.A.  Anderson,  G.  Moddel,  M.A.  Paesler,  and  H.  Paul,  "The  Importance 
of  Argon  Pressure  in  the  Preparation  of  r.f.  Sputtered  Amorphous  Si-H 
Alloys,"  to  be  published  in  J.  Vac.  Science  and  Technology. 


II. 11  Vacancy  Stabilization.  L.  Huisman,  A.  Carlsson,  H.  Ehrenreich,  and 
C.D.  Gelatt,  Contracts  N0014-75-C-0648,  NSF  DMR77-24295  and  NSF 
DMR77-10210;  Research  Unit  1. 

First  principles  calculations  of  the  effect  of  vacancies  on  the  heat 
of  formation  of  the  refractory  compounds  TiO  and  TiC  have  provided  insight 
into  the  qualitative  features  which  distinguish  the  vacancy-forming 
properties  of  these  systems.1  In  stoichiometric  TiC  at  low  temperature, 
the  equilibrium  vacancy  concentration  is  less  than  one  percent,  while  in 
TiO  it  is  on  the  order  of  ten  percent  on  both  sublattices.  Detailed 
calculations  of  the  electronic  structure,  and  simple  model  Hamiltonian 
analogs,  demonstrate  the  features  of  the  metal -vacancy  interaction  which 
lead  to  vacancy  stabilization.  This  new  insight  should  prove  useful  in 
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understanding  defect  formation  in  other  systems,  such  as  reconstructed 
surfaces  and  distorted  phases  of  nominally  A-15  structure  compounds,  that 
exhibit  defects  in  a high-symmetry  structure  at  equilibrium. 


Reference 


1.  L.  Huisman,  A.  Carlsson,  H.  Ehrenreich,  and  C.D.  Gelatt,  "Defects 

Systems:  A New  Stabilization  Mechanism  in  Compounds,"  to  be  published. 
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11.12  Electronic  Theory  of  Disordered  Semiconductors.  C.D.  Gelatt  and 

S.  Kivelson,  Contracts  N00014-75-C-0648  and  NSF  DMR77-24295; 

Research  Unit  1. 

The  report  of  the  effective  mass  theory  formalism  for  non-crystalline 

systems  is  completed.*  We  find  that  charged  substitutional  impurities  in 

a model  of  amorphous  Si  induce  large  radius,  small  binding  energy  impurity 

states  whose  properties  are  quite  similar  to  those  of  the  corresponding 

impurity  states  in  crystalline  Si.  The  qualitative  features  of  this 

2 

calculation  have  been  confirmed  in  exactly  soluble  extension  of  the  Lloyd 
model.  The  influence  of  such  spatially  large  states  on  the  temperature  and 
time  dependence  of  photoluminescence  in  amorphous  Si  is  currently  under 
investigation . 

We  have  analyzed  conductivity  in  disordered  insulators  in  which 
carriers  move  by  phonon-assisted  hopping,  and  demonstrated  that  the  DC 
conductivity  is  finite,  even  in  systems  characterized  by  a large  dispersion 
in  hopping  integrals  among  sites  with  a narrow  range  of  energies,  in 
contrast  to  the  vanishing  DC  conductivity  predicted  by  the  results  of 
Scher  and  Montroll  based  on  a single-site  decoupling  of  the  configuration 
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average  developed  by  Scher  and  Lax.  This  result  calls  into  question 
conclusions  about  the  electronic  states  in  such  systems  drawn  from  analysis 
of  transient  photoconductivity  using  the  Scher  and  Montroll  formalism. 


References 


1.  S.  Kivelson  and  C.D.  Gelatt,  Jr.,  "Effective  Mass  Theory  in  Non- 
Crystalline  Solids,"  Phy8.  Rev.  B,  to  be  published. 

2.  S.  Kivelson  and  C.D.  Gelatt,  Jr.,  "Impurity  States  in  a Disordered 
Insulator:  The  Lloyd  Model,"  submitted  to  Pkya.  Rev.  B. 

3.  S.  Kivelson,  "Hopping  Transport  and  the  Continuous-Time  Random  Walk," 
submitted  to  Phys.  Rev.  B. 


11.13  Interfaces . A.  Carlsson  and  H.  Ehrenreich,  Contract  N00014-75-C- 
0648;  Research  Unit  1. 


Ehrenreich' s interests  in  semiconductor  theory  have  been  renewed 
because  of  his  involvement  as  Chairman  of  a group  formed  at  the  suggestion 
of  the  White  House  OSTP  studying  solar  photovoltaic  energy  conversion.* 

As  a member  of  ARPA's  Materials  Research  Council  he  has  been  and  continues 
to  be  involved  in  investigations  concerning  physical  size  limitations  in 
computer  elements  in  connection  with  very  large  scale  integration. 

A better  understanding  of  interfaces  between  semiconductors  and 
oxides,  semiconductors  and  metal,  and  even  semiconductors  and  electrolytes 
represents  an  important  ingredient  in  both  areas.  Questions  such  as  the 
conditions  under  which  one  expects  states  localized  in  the  vicinity  of  the 
interface  and  the  effects  of  disorder  on  interface  states  are  being  studied 
by  the  use  of  model  Hamiltonians . The  Heine-Weaire  Hamiltonian  has  been 
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used  to  study  the  electronic  structure  of  the  interface  between  a semi- 
conductor and  simple  cubium  metal.  This  same  Hamiltonian  has  also  been 
generalized  to  apply  to  SiO^  in  order  to  investigate  some  of  the  basic 
properties  of  the  Si-Si02  interface,  in  particular,  the  formation  of 
localized  states.  It  is  hoped  that  a better  understanding  of  the 
passivation  process  will  emerge  from  the  latter  investigation  and  that 
additional  insight  will  be  obtained  concerning  the  height  of  Schottky 
barriers  from  studies  of  the  semiconductor-metal  interface. 

Reference 

1.  Principal  Conclusions  of  the  American  Physical  Society  Study  Group  on 
Solar  Photovoltaic  Energy  Conversion  (H.  Ehrenreich,  chairman) , 
American  Physical  Society,  1979. 


11.14  Electrons  in  a Random  Potential.  M.J.  Stephen,  Contract  N00014-75- 

C-0648?  Research  Unit  1. 

According  to  the  band  theory  an  electron  moving  in  a random  potential 
is  described  by  a diffusion  equation.  However,  for  sufficiently  large 
disorder,  diffusion  becomes  impossible  and  the  electron  becomes  localized 
in  space.  A critical  energy  Ec  (the  mobility  edge)  in  the  band  of  allowed 
states  separates  localized  states  from  extended  states  and,  depending  on 
the  position  of  the  Fermi  energy,  the  material  behaves  as  a semiconductor 
or  a metal.  We  have  developed  a self-consistent  theory  to  describe  the 
mobility  edge.  Self-consistent  integral  equations  for  the  distribution 
functions  of  the  real  and  imaginary  parts  of  the  electron  self-energy  are 
obtained.  The  self-consistent  approximation  consists  in  neglecting  the 
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effects  of  closed  loops.  The  theory  is  exact  in  one  dimension  and  is 
equivalent  to  a mean  field  approximation  in  higher  dimensions.  The 
determination  of  the  mobility  edge  reduces  to  the  solution  of  a linear 
integral  equation  and  this  has  been  solved  approximately  for  some  choices 
of  the  random  potential.  The  determination  of  the  diffusion  constant  and 
i correlation  length  requires  a solution  of  the  integral  equation  near  the 

mobility  edge  and  this  work  is  in  progress. 
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III.l  Data  Base  System  Problems.  Y.C.  Ho,  B.  Kurtaran,  and  D.M.  Chiu, 
Contracts  N00014-75-C-0648,  N00014-78-C-0531 , and  HSF  Grant 
ENG  76-11824;  Research  Unit  9. 

(i)  Query  Optimization 

Consider  a single  query  such  as  "how  many  foreign  born  males  over  35 
are  teaching  in  an  eastern  university  with  a student  population  over  15,000 

I 

and  whose  research  interest  is  in  system  theory?"  To  retrieve  the  answer 
from  the  U.S.  census  data  base,  one  can  visualize  a pass  through  the  census 


data  with  an  appropriate  filter.  However,  if  one  is  to  answer  hundreds  and 
thousands  of  such  queries  per  day  using  the  data  base,  then  such  a naive 


approach  is  clearly  infeasible.  Optimization  in  terms  of  file  management, 
storage  allocation,  and  the  sequencing  and  processing  of  intermediate  results 

become  necessary.  Two  pertinent  and  recent  references  are  Hong  and 

1 2 
Youssefi  and  Eswaran . 

Abstractly,  the  problem  can  be  considered  as  follows:  Let  the  vector 
x represent  the  "state"  of  the  data  base.  It  is  supposed  to  reflect  the 
size  of  the  various  relations  currently  in  store,  their  access  mechanisms, 
indices  and  directories,  and  all  the  logical/physical  parameters  effective 
retrieval.  As  we  process  the  query,  intermediate  results  and  relations 
appear  which  take  the  data  base  to  a new  "state"  which  is  determined  by  the 
choice  of  intermediate  processing  decisions.  Thus,  a dynamic  control  problem 
takes  shape.  The  object  is  to  choose  the  processing  decision  optimally  so 
that  some  criterion  such  as  the  processing  time  or  effort  is  minimized 
subject  to  the  constraints  that  the  final  state  contains  the  desired  output. 
We  have  formulated  two  simplified  versions  of  the  query  optimization  problem 
in  a distributed  data  base.  A report  is  under  preparation. 


(ii)  Security  and  Access 

Typically,  in  data  base  systems,  the  concept  of  an  access  matrix 
(rows  represent  users;  columns,  data)  is  used  to  determine  "who  can  access 
what".  As  various  "grant"  and  "revoke"  coimnands  are  given  by  different 
users  to  each  other  in  a distributed  and  random  manner,  the  access  matrix 
or  the  state  of  access  evolves  in  a rather  complex  way  effecting  the 
security  levels  (exposure)  or  various  data  items.  What  is  generally  desired 
is  that  some  overall  security  requirement,  such  as  "no  low  level  user  should 
gain  access  to  high  level  data"  be  insured.  To  use  a gross  analogy,  by 
insisting  that  Re  A. (A)  <0,  a local  property,  we  can  insure  that  x(t) 
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(defined  through  x * Ax)  approaches  zero  as  t-*-00,  a global  property. 
Similarly,  by  imposing  local  requirements  on  how  access  is  modified  by  local 
transition  of  "grant"  and  "revoke",  we  wish  to  insure  overall  secure 
information  flows.  Two  pertinent  references  are  (3)  and  (4).  Our  work5 
represents  a marriage  of  the  ideas  in  (3)  and  (4). 

J 
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III. 2 Teams  and  Organizations.  Y.C.  Ho  and  N.  Papadopoulos , Contracts 

N00014-7 5-C-0648 , N00014-78-C-0531 , N00014-77-C-0533,  and  NSF  Grant 
ENG  76-11824;  Research  Unit  9. 

We  studied  the  question  of  "who  should  know  what"  in  the  context  of  a 
general  quadratic  team  problem  with  linear  information  structure.  In  this 
context  two  matrices  C and  H+  are  of  central  importance.  The  matrix  C 
is  dependent  only  on  the  payoff  structure  while  H+,  being  the  generalized 


inverse  of  the  information  structure  matrix  H,  depends  only  on  it.  Rows  of 


C shows  how  different  states  of  nature  (i.e.,  the  environment)  effects 
each  decision  maker.  Columns  of  H+  detail  how  each  piece  of  information 
can  help  in  estimating  the  states  of  nature.  The  combination  of  C and  H 
then  answers  the  question  "which  piece  of  information  is  of  use  to  which 
decision  maker".  The  implication  of  this  to  the  organisation  of 
communication  is  obvious. 

Reference 

1.  Y.C.  Ho  and  N.  Papadopoulos , "Further  Notes  of  Redundancy  in  Teams," 
IEEE  Transactions  on  Automatic  Control,  April  1979  (to  appear) . 


III. 3 Buffer  Storage  Problem  in  a Production  Line.  Y.C.  Ho  and  A.  Eyler, 
Contract  N00014-75-C-0648,  and  The  Charles  Stark  Draper  Lab.,  Inc., 
Contract  DL-H-1 62176;  Research  Unit  9. 

A problem  of  long  standing  in  manufacturing  is  the  optimal  allocation 
of  buffer  storage  in  a serial  production  line  which  is  subject  to  random 
failures  and  repairs  of  its  component  machines.  This  problem  can  be 
formulated  as  a problem  in  Markov  chain  analysis  with  combinatorially  large 
state  space.  We  present  a general  solution  to  this  problem.  The  key 
ingredient  in  the  solution  is  an  efficient  scheme  for  the  computation  of  the 
gradient  component  of  the  line  throughput  to  buffer  size  increases  at 
individual  locations.  Analytical  results  and  experimental  verification 
have  been  obtained.  The  method  of  gradient  calculation  is  general  and  can 
be  adopted  to  calculate  sensitivities  of  production  efficiency  for  other 


parameters . 
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1.  Y.C.  Ho,  A.  Eyler,  and  T.T.  Chien,  "A  Gradient  Technique  for  General 
Buffer  Storage  Design  in  a Serial  Production  Line,”  Prc^^edinge  of  the 
1978  IEEE  Decision  and  Control  Conference , San  Diego,  Calif.,  January 
1979. 


III. 4 Decentralized  Resource  Management  in  Large  Systems.  Y.C.  Ho  and 

R.  Suri,  Contract  N00014-75-C-0648,  NSF  Grant  ENG  76-11824,  and  The 
Charles  Stark  Draper  Lab.,  Inc.,  Contract  DL-H-162176;  Research  Unit  9. 

The  study  of  Large-Scale  Systems  is  becoming  increasingly  important, 

as  is  reflected  by  the  recent  devotion  of  an  entire  issue  of  a prominent 

journal1  to  the  topic  of  "Large-Scale  Systems  and  Decentralized  Control." 

2 3 

Our  recent  work  ' studies  the  application  of  decentralization  to  one  aspect 
of  such  systems,  namely , the  problem  of  Resource  Management  in  Large  Systems. 

In  a large  operational  system,  where  a very  large  number  of  activities 
share  a number  of  limited  resources,  this  Resource  Management  problem  has 
three  main  objectives.  The  first  (the  "Initial  Allocation"  or  "Design" 
problem)  is  to  find  an  assignment  of  resources  to  all  the  activities,  such  that 
all  the  system  constraints  are  satisfied,  and  all  activities  are  operating,  so 
that  the  system  as  a whole  can  operate'.  The  second  (the  "New-Assignment" 
problem)  is  to  find  a rationale  for  allocating  resources  to  new  activities. 

It  is  presumed  that  new  activities  are  initiated  frequently  enough  that  we 
do  not  wish  to  resolve  the  entire  problem.  The  third  objective  ("Periodic 
Review"  problem)  is  to  find  an  efficient  way  of  reallocating  resources  in 
order  to  reflect  the  changing  needs  of  the  individual  activities,  as  well  as 
the  changes  in  total  resource  usages. 


II 1-6 


The  main  contributions  of  our  work  is  to  extend  the  applicability 
of  decentralized  techniques  to  problems  where  the  resource  usages  functions 
are  ill-behaved  (e.g.,  discontinuous).  We  prove  the  existence  of  a 
decentralized  solution,  and  give  efficient  algorithms  guaranteed  to  converge 
to  this  solution. 

Our  approach  enables  efficient  solution  of  a problem  with  over 
60,000  activities,  each  of  which  has  a 16-dimensional  decision  vector. 

This  one-million  dimensional  problem  would  be  hard  and  impractical  to  solve 
by  other  methods. 
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III. 5 Identification  of  Time  Invariant  Linear  Systems  from  Non-Station 


Cross-Sectional  Data.  P.E.  Caines,  Contract  N00014-75-C-0648; 
Research  Unit  10. 


We  have  developed  a new  theory  and  methodology  for  the  identification 
of  linear  time  invariant  finite  dimensional  systems  from  non-stationary 
cross-sectional  data.  The  parameters  of  a system  model  are  estimated  from 
repeated  observations  on  a set  of  transient  output  samples  of  the  system, 
where  each  sample  is  of  finite  duration.  This  work  and  its  underlying 
theory  is  described  in  the  papers  (1-4) ; it  permits  parameter  identification 
for  a class  of  systems  for  which  this  was  not  previously  possible. 
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In  our  theory  we  formulate  the  notion  of  (non-stationary)  system 

. 2 

identif iability  as  that  of  unique  stochastic  realizability.  That  is  to  say, 

a parameterized  system  is  identifiable  if  (exact)  estimates  of  the  output 

process  covariance  lead  to  unique  estimated  values  of  the  system  parameters. 

2 

We  have  recently  shown  that  for  non-stationary  systems  of  the  form 
appearing  in  (1)  one  has  (local  identifiability)  without  any  condition  on 
the  stability  of  the  system.  This  result  gives  significant  scope  to  the 
applications  of  the  technique  described  in  (1) . In  addition  the  computer 
program  presented  in  (5)  is  being  currently  improved  in  terms  of  its 
efficiency  and  flexibility. 
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III. 6 Deterministic  and  Stochastic  Adaptive  Control.  P.E.  Caines, 


Contract  N00014-75-C-0648;  Research  Unit  10. 


We  have  studied  the  problem  known  in  systems  science  as  the  adaptive 
control  problem.  The  issue  here  is  whether  it  is  possible  to  stabilize  an 
unstable  (linear)  dynamical  system  whose  parameters  are  unknown.  This  is 
to  be  achieved  by  simultaneously  estimating  the  system  parameters  and 
continuously  designing  improved  feedback  control  laws.  Although  there 
exist  several  examples  of  successfully  operating  adaptive  control  regulators 
in  industry,  a theoretical  study  of  their  behavior  has  eluded  research  workers 
for  years.  Some  important  partial  results  have  recently  appeared  in  (1)  and 


We  showed  that  stabilizing  adaptive  control  is  indeed  possible  for 
multi-input  multi-output  linear  systems  with  the  control  laws  being  chosen 
from  a large  class  of  simple  and  intuitively  appealing  parameter  estimation 
and  control  algorithms.  In  addition,  by  the  use  of  recent  results  in 
probability  theory,  the  analogous  problem  for  stochastic  systems  was  solved. 
It  should  be  emphasized  that  although  the  underlying  systems  to  be  controlled 
are  linear  the  closed  loop  system  with  its  learning  and  control  regulator  is 
nonlinear. 
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III. 7 Nonlinear  Systems.  R.W.  Brockett,  Contract  N00014-75-C-0648;  Research 

Unit  8. 

In  our  paper1  we  examine  the  structure  of  the  necessary  conditions  of 
optimal  control  using  the  Volterra  series  results  developed  over  the  last 
few  years.  This  approach  goes  beyond  earlier  efforts  in  this  direction  in 
that  it  includes  the  effect  of  terminal  constraints  and  describes  the  minimal 
amount  of  smoothness  on  the  system  and  performance  measure  necessary  to  find 
necessary  conditions  which  are  also  sufficient.  The  results  of  this  paper 
are  paralleled  by  recent  developments  in  the  mathematical  programming 
literature  associated  with  the  group  around  Ben-Israel. 
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III. 8 Estimation  Theory.  R.W.  Brockett,  Contract  N00014-75-C-0648; 
Research  Unit  8. 


In  reference  (1)  we  describe  a new  point  of  view  on  the  conditioned, 
density  equation — the  basic  equation  of  nonlinear  filtering.  We  show  there 
the  connection  between  the  existence  of  simple  filters  and  certain  Lie 
algebras.  This  point  of  view  allows  us  to  describe  certain  estimation 
problems  which  are  not  linear  and  Gaussian  but  which  have,  nonetheless, 
simple  filtering  algorithms. 
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III. 9 Linear  Systems.  R.W.  Brockett,  Contract  N00014-75-C-0648;  Research 
Unit  8. 


We  have  just  completed  work  on  the  celebrated  "output-feedback" 
problem,  obtaining  by  classical  geometric  techniques  the  most  extensive, 
and  in  some  cases  the  sharpest,  results  currently  available.  This  has 
ramifications  in  stability  theory  and  the  same  techniques  also  apply  to 
problems  of  adaptive  control,  e.g.,  the  stabilization  of  systems  depending 
on  parameters.  In  reference  (1)  we  discuss  entirely  unexpected  connection 
between  the  theory  of  linear  and  nonlinear  systems. 
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Research  in  the  area  of  electromagnetic  radiation  is  directed  toward 
the  solution  of  practical  problems  through  the  complete  understanding  of 
the  underlying  physical  phenomena.  This  involves  the  coordinated  application 
of  modern  analytical,  numerical,  and  experimental  techniques  and  the  use  of 
high-speed  computers  and  precision  instrumentation.  Most  practically 
significant  problems  in  this  area  are  sufficiently  complicated  that  extensive 
computation  and  measurement  are  often  required  to  justify  approximations  that 
are  usually  necessary.  Where  possible,  general  formulas  are  obtained  and 
verified  experimentally  so  that  the  phenomenon  under  study  can  be  understood 
physically  in  analytical  form  and  not  just  as  a set  of  numbers. 

The  researches  are  conveniently  grouped  under  two  main  headings,  viz., 
cmtennaa  and  waves  in  and  over  dielectric t conducting , and  plasma  media, 
and  antennas  in  air.  The  first  group  is  concerned  with  the  circuit, 
radiating,  and  scattering  properties  of  antennas  and  arrays  and  the  waves 


generated  by  them  in  their  dependence  on  the  electrical  properties  o£  a 

1 ' I 

material  medium  (lake  water,  sea  water,  earth,  ice,  living  tissue,  plasma, 

etc.)  in  various  forms  and  degrees  of  proximity.  The  antennas  may  be 

* 

directly  embedded  in  the  medium  or  separated  from  it  by  a layer  of  in-  f 

if  * 

sulation.  The  second  group  of  investigations  deals  with  the  detailed 

* 
b 

electrical  properties  of  complicated  metal  radiating  and  scattering 
structures  in  air  including  especially  electrically  thin  and  thick  crossed 
conductors  used  to  model  an  aircraft. 
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ANTENNAS  AND  WAVES  IN  AND  OVER  DIELECTRIC, 


CONDUCTING,  AND  PLASMA  MEDIA 


IV. 1 Theoretical  and  Experimental  Study  of  the  Scattering  from  an  Obstacle 

Above  the  Earth.  H.-M.  Lee,  T.T.  Hu,  and  L.C.  Shen,  Contracts 

I 

N00014-7 5-C-0648  and  F19628-77-C-0146;  Research  Unit  11. 

A knowledge  of  the  electromagnetic  field  scattered  by  targets  located 
on  or  above  the  surface  of  the  earth  is  of  primary  importance  in  radar. 

This  field  is  radiated  by  the  currents  induced  in  the  target  by  the  incident 
plane  wave.  It  follows  that  a first  step  in  evaluating  the  scattered  field 
is  the  determination  of  the  induced  currents.  These  depend  not  only  on  the 
geometrical  and  electrical  characteristics  of  the  scattering  target  and  the 
direction  and  polarization  of  the  incident  wave,  but  also  on  the  proximity 
of  the  earth  with  its  conducting  and  dielectric  properties.  In  effect,  the 
current  in  the  target  corresponds  to  a continuous  distribution  of  infinites- 
imal dipoles,  each  with  appropriate  amplitude  and  relative  phase.  Thus,  the 
toted,  scattered  field  is  the  vector  sum  of  the  contributions  by  all  of  the 
small  dipoles,  each  with  its  particular  amplitude  and  phase.  Hie  field  of 
one  of  the  infinitesimal  dipoles  is  the  solution  of  the  Sommerfeld  problem. 

The  thin  circular  loop  is  the  most  convenient  target  to  analyze 

1 j 

initially  because  of  its  rotational  symmetry  and  its  usefulness  in  the 
construction  of  composite  targets.  Furthermore,  the  current  distribution 
and  scattered  field  are  well  known  when  the  loop  is  isolated,  i.e., 
sufficiently  far  from  the  earth.  Among  loops  that  possess  rotational 
symmetry  around  the  z-axis  normal  to  the  earth's  surface,  the  cylindrical 
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loop  (cut  from  a cylinder  of  zero  thickness  with  its  height  much  smaller 

than  its  radius  and  the  incident  wavelength)  has  the  simplest  integral 

equations  for  the  two  components  of  the  surface  current  excited  by  an 

external  field.  Although  the  two  integral  equations  can  be  decoupled, 

their  evaluation  is  exceedingly  difficult  due  to  the  many  widely  varying 

parameters . In  a recent  report,3  a variational  approach  is  used  to  solve 

the  problem.  The  solution  is  extended  to  multiple  loops  and  to  multiple 

loops  each  with  arbitrary  cross  section.  Attempts  are  continuing  to  obtain 

an  analytical  solution  to  the  decoupled  integral  equations.  The  lowest- 

order  analytical  solution  to  this  problem  was  originally  derived  by  T.T.  Hu 

2 

in  his  study  of  the  loop  antenna.  The  next -order  correction  to  the  kernel 
of  the  decoupled  integral  equations  has  now  been  obtained.  It  is  expected 
that  this  knowledge  will  help  in  the  solution  of  other  problems,  including 
the  loop  placed  above  the  earth. 

The  circular  loop  is  also  a simple  target  for  use  in  a systematic 

program  of  measurement.  As  described  in  detail  in  am  earlier  report,3  one 

set  of  measurements  has  been  completed  of  the  field  scattered  by  horizontal 

resonant  and  nonresonant  circular  loops  and  by  a straight  wire  one-half 

wavelength  in  length  over  a wide  range  of  sizes,  angles  of  incidence,  and 

heights  outdoors  over  the  earth.  These  results  were  also  presented  as  a 
4 

conference  paper.  An  apparatus  has  been  constructed  to  continue  the 
measurements  systematically  over  a tank  of  water  with  variable  conductivity. 

A comprehensive  experimental  study  is  in  progress. 
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IV. 2 Numerical  and  Analytical  Determination  of  Electromagnetic  Fields  of 
Antennas  in  the  Earth  Near  Its  Interface  with  Air.  R.W.P.  King, 

B.H.  Sandler,  and  L.C.  Shen,  Contract  N00014-75-C-0648;  also  NSF 
Grant  ENG-76-17596  with  the  University  of  Houston;  Research  Unit  11. 

The  calculation  of  the  electromagnetic  field  in  a dissipative  half- 
space due  to  an  infinitesimal  electric  dipole  in  the  same  half -space  and 
with  its  axis  parallel  to  the  interface  is  a classic  problem  closely  related 
to  the  work  of  Sommerfeld  on  the  field  in  air  due  to  a dipole  also  in  air 
above  a dissipative  earth.  It  has  been  the  subject  of  extensive  investi- 
gations. In  an  earlier  paper,1  expressions  for  the  three  cylindrical 
components  of  the  electric  field  were  obtained  from  Bafios*  general  formulas 
in  terms  of  complex  integrals.  These  were  evaluated  numerically  for  a 
range  of  frequencies  and  radial  distances  from  a horizontal  dipole  in  sea 

water,  lake  water,  and  dry  earth.  An  alternative  derivation  that  proceeds 

2 

directly  from  Maxwell's  equations  and  the  boundary  conditions  has  led  to 
equivalent  integral  expressions.  These  have  been  evaluated  by  a less  time- 
consuming  program.  As  a result,  more  comprehensive  calculations  have  been 
possible.  Tables  have  been  prepared  of  the  generally  useful,  radial 
component  of  the  electric  field  for  radial  distances  from  100  m to  100  km. 


*¥* 


properties  of  the  earth.  The  field  generated  by  the  currents  in  an  extended 

conductor  can  be  calculated  from  the  known  field  of  am  infinitesimal  dipole 

only  after  the  distributions  of  current  have  been  determined  in  their 

4 

dependence  on  the  properties  of  the  earth.  A paper  has  been  written  which 
treats  three  special  cases  of  a horizontal  wire  in  or  close  to  a dissipative 
earth.  These  are  an  end-driven  Beverage-type  antenna,  a horizontal  wire  as 
a scatterer  in  the  presence  of  the  earth,  and  an  embedded  insulated 
travel ing-wave  antenna.  B'or  each  of  these,  an  analytically  determined 
distribution  of  current  is  available  that  takes  full  account  of  the 
electrical  properties  of  the  adjacent  or  surrounding  material  half-space 
provided  the  medium  is  electrically  dense  and  the  height  of  the  antenna 
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(when  in  air)  is  sufficiently  small.  For  the  embedded  insulated  antenna 

with  insulation  radius  b,  the  depth  d below  the  surface  must  be  large 

2 2 

enough  so  that  the  medium  is  effectively  infinite,  i.e.,  d »b  . Also, 
the  insulated  part  of  the  antenna  must  be  electrically  thin  in  terms  of 
both  wave  numbers  of  the  material  half-space  and  kQ  of  the  insulation, 

viz.,  kjb<l  and  kQb  < 1.  For  each  of  the  three  cases,  expressions  for 
the  most  useful  component  of  the  radiated  or  scattered  field  are  formulated. 

In  a related  study6  special  attention  is  given  to  the  eccentrically 
insulated  traveling-wave  antenna  embedded  beneath  the  surface  of  the  earth. 
This  antenna  is  particularly  useful  for  subsurface  communication  by  lateral 
waves  since  its  field  pattern  is  highly  directive  upward.  Single  analytical 
formulas  for  the  input  impedance,  current  distribution,  and  field  pattern 
are  derived  that  incorporate  the  recent  advances  in  determining  the 
properties  of  antennas  of  finite  size  in  the  presence  of  the  air-earth 
boundary.  A step-by-step  procedure  for  designing  an  eccentrically  insulated 
traveling-wave  antenna  is  given.  To  increase  the  directivity  further,  a 
broadside  array  of  eccentrically  insulated  traveling-wave  antennas  can  be 
used. 

Since  lateral-wave  transmission  is  inefficient  even  when  fairly 
directive  eccentrically  insulated  traveling-wave  antennas  are  used, 
relatively  large  amounts  of  power  must  be  supplied  to  the  antenna  in  order 
to  generate  a significant  field  at  a distance.  The  question  arises:  How 
intense  is  the  field  radiated  upward  into  the  air?  Is  this  sufficiently 
intense  to  be  hazardous  to  living  organisms?  This  is  the  subject  of  a new 
paper,6  in  which  accurate  calculations  cure  made  of  the  electric  and  magnetic 
fields  in  the  air  above  a horizontal  electric  dipole  in  the  sea,  a lake,  and 
the  earth.  This  is  carried  out  directly  from  the  exact  general  integrals. 
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Useful  approximate  expressions  for  the  field  in  air  above  the  transmitting 
dipole  are  also  obtained  which  indicate  that  the  limiting  forms  given  by 
Bafios  are  oversimplified  for  the  determination  of  the  phases.  Numerical 
calculations  from  the  exact  formulas  are  also  made  of  the  electric  and 
magnetic  fields  on  the  surface  directly  over  the  embedded  dipole,  and 
approximate  simple  formulas  cure  derived  and  compared.  In  all  cases,  the 
Poynting  vector  is  directed  upward  into  the  air  and  not,  as  stated  by  Bafios, 
downward  from  the  air  into  the  conducting  half-space.  These  results  cure  in 
complete  agreement  with  the  theory  of  embedded  antennas  which  indicates  that 
it  is  possible  to  transmit  from  a dipole  within  a dielectric  or  conducting 
body  out  into  the  surrounding  air,  - a process  which  requires  an  upward 
directed  Poynting  vector  in  the  air. 
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IV. 3 Transmission  and  Reception  with  Bare  and  Insulated  Antennas  Embedded 

in  Skin.  R.W.P.  King  and  B.H.  Sandler,  Contract  N00014-75-C-0648;  and  NSF 

Grant  ENG  73-03507  A01  with  Northeastern  University;  Research  Unit  11. 

Transponders  are  embedded  in  or  located  just  outside  a living  organism 
for  the  remote  monitoring  of  various  quantities  such  as  body  temperature, 
heart  beat,  and  intensity  of  radiation.  Alternatively,  it  may  be  of 
interest  to  use  an  embedded  or  externally  attached  transponder  to  identify 
or  locate  animals  in  their  environment.  The  successful  operation  of  such 
systems  depends  greatly  on  the  proper  design  of  transmitting  and  receiving 
antennas  for  use  within  or  near  a living  organism.  For  this,  it  is 
necessary  to  know  the  magnitude,  phase  and  polarization  of  the  electric 
field  at  different  points  of  observation  in  the  body  and  just  outside  its 
surface.  Furthermore,  the  complete  circuit  properties  (admittance,  current, 
and  charge  per  unit  length)  of  the  embedded  or  attached  antenna  are  of 
interest  in  their  dependences  on  the  electrical  properties  of  the  environment. 

The  electric  field  in  a body  has  the  commonly  assumed  simple  properties 
of  an  exponentially  attenuated  traveling  plane  wave  only  when  (a)  the 
exciting  field  is  a normally  incident,  plane  wave,  (b)  the  body  is  large 
enough  to  permit  its  surface  to  be  treated  as  locally  plane,  and  (c)  internal 
reflections  can  be  neglected.  At  other  than  normal  incidence,  tr.e  wave 
transmitted  into  the  body  is  inhomogeneous  with  surfaces  of  constant  phase 
and  constant  amplitude  that  neither  coincide  nor  travel  in  the  same  direction. 

It  follows  that  the  field  acting  to  induce  currents  in  an  embedded  antenna 
is  quite  different  from  that  acting  on  the  same  antenna  in  air  at  the  same 
other  than  normal  angle  of  incidence.  The  determination  of  the  current  is 
further  complicated  at  any  angle  of  incidence  if  the  attenuation  constant 
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of  the  medium  is  low  enough  or  if  the  body  is  small  enough  or  consists  of 
electrically  different  layers  to  make  internal  reflections  and  the 
associated  standing  waves  significant . 

A complete  investigation  has  been  reported3  on  the  voltage  generated 

r 

by  an  incident  plane  wave  across  the  load  of  a bare  or  insulated  receiving 
dipole  at  various  depths  in  a half-space  composed  of  a single  homogeneous  * 

material  with  the  electrical  properties  of  skin.  This  includes  the  deter- 
mination of  the  impedance  of  the  antenna  in  the  particular  environment. 

The  reciprocal  problem  of  obtaining  the  voltage  across  the  load  in  a 
receiving  antenna  located  at  a distant  point  in  air  by  an  emf  applied  to  a 

bare  or  insulated  dipole  embedded  at  various  depths  in  a half-space  of  skin 

2 

has  also  been  solved.  Specifically,  two-way  communication  between  an 
element  embedded  up  to  1 cm  in  skin  and  an  outside  unit  at  a distance  of 
4.3  m at  the  optimum  frequency  of  700  MHz  or  1 m at  3 GHz  has  been  analyzed. 

The  solution  of  the  electromagnetic  problem  of  the  transponder  at  various 
depths  in  a homogeneous  half-space  of  skin  is  thus  completed. 

A further  investigation  has  been  made3  to  extend  these  analyses  to 
receiving  and  transmitting  antennas  located  in  or  near  the  surface  of  a 
three-layered  half-spacing  consisting  of  finite  thicknesses  of  skin  and 
fat  on  an  infinite  depth  of  muscle.  In  particular,  the  properties  of  bare 
and  insulated  antennas  are  examined  when  successively  located  in  the  air 
just  outside  the  body  and  at  the  centers  of  the  layers  of  skin  and  fat, 
and  at  a fixed  distance  in  the  muscle.  The  electromagnetic  problems  that 
require  solution  for  the  centrally  loaded  receiving  antenna  include  the 
determination  of  the  following  quantities:  (a)  The  electric  field  at  all 
points  where  a transponder  may  be  located  in  terms  of  the  incident  field 
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at  the  surface  of  the  skin,  (b)  The  driving-point  impedance  ZQ  of  the 
bare  and  insulated  antennas  when  located  successively  in  the  air  and  in 
the  layers  of  skin,  fat  and  muscle.  (c)  The  current  in  the  load  for  both 
types  of  antenna  in  each  of  the  four  locations.  (d)  The  voltage  V 

L 

across  the  load  Z . For  the  transmitting  mode,  the  additional  quantity 
required  is  the  ratio  |vl/|ve|  of  the  voltage  V.  across  the  load  of 

Ii  Li 

0 

the  receiving  antenna  in  air  to  the  driving  voltage  V applied  across 
the  terminals  of  the  embedded  transmitting  antenna.  Bare  and  insulated 
antennas  are  considered  over  a frequency  range  up  to  3 GHz  with  layer 
thicknesses  of  skin  from  0.1  to  1.0  cm,  and  of  fat  from  0 to  1.5  cm.  The 
application  of  the  results  to  finite  bodies  is  discussed. 
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IV. 4 Theoretical  ar.d  Experimental  Study  of  Electromagnetic  Fields  and 
Antennas  in  Dissipative  and  Dielectric  Cylinders.  R.  Bansal, 

Contract  N00014-75-C-0648;  Research  Unit  11. 

Tne  complex  physical  problem  associated  with  the  transmission  of  an 
electromagnetic  signal  from  an  antenna  in  air  to  another  located  at  an 
arbitrary  distance  inside  a finite  material  body  is  of  considerable  practical 
interest.  Such  an  investigation  is  of  importance,  for  example,  in  the  study 
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of  the  properties  of  antennas  embedded  in  dielectric,  poorly  conducting 
aircraft,  in  the  determination  of  the  biomedical  hazards  of  EM  radiation 
of  living  biological  systems,  and  in  the  use  of  embedded  transponders. 

The  problem  is  also  closely  related  to  diffraction  by  bodies  in  the  resonant 

r | 

regions  and  to  scattering  by  dielectric  objects.  * 

As  a first  step  in  the  investigation,  a systematic  experimental  study  * 

has  been  initiated  to  determine  the  distribution  of  the  electric  field  (in 
both  amplitude  and  phase)  inside  a finite  circular  cylinder  of  water 
enclosed  in  a styrofoam  container  when  illuminated  by  an  approximately  plane 
wave.  An  experimental  setup  for  this  purpose  has  been  designed,  fabricated, 
and  thoroughly  tested.  The  setup  includes  an  indoor  aluminum  ground  plane, 
a transmitter  consisting  of  a quarter-wave  monopole  backed  by  a corner- 
reflector,  and  the  probing  system  located  at  the  receiving  end  inside  a 
water  cylinder.  For  experimental  reasons  a thin  central  conducting  tube  is 
included  along  the  axis  of  the  cylinder  of  water;  this  is  needed  to  shield 
the  transmission  lines  leading  to  the  probes.  The  probing  system  makes  use 
of  a pair  of  dipoles  to  measure  simultaneously  the  amplitude  and  phase  of 
the  axial  and  transverse  components  (the  z and  <j>  components  in  the 
usual  cylindrical  coordinates)  of  the  electric  field.  The  radial  component 
(p  component)  of  the  electric  field  cannot  be  measured  in  the  present  setup 
since  the  radial  arm  supporting  the  dipoles  necessarily  interferes  with  the 
radial  component  of  the  electric  field.  r 

At  low  frequencies  it  is  quite  difficult  to  generate  an  incident 

J 

4ft  •' 

plane-wave  field  indoors.  On  the  other  hand,  since  the  wavelength  in  water 
is  much  shorter  than  the  free-space  wavelength,  it  is  undesirable — from 
physical  design  considerations  in  the  receiving  system — to  use  very  high 
transmission  frequencies.  As  a compromise  between  these  two  conflicting 

i I 

— ~ . ..  r .j-.-iji 


IV-13 


requirements,  10  Hertz  has  been  selected  as  the  transmitter  frequency,  at 
least  initially.  The  setup  can  be  modified  easily  to  make  measurements 
over  a wide  band  of  frequencies. 

After  the  preliminary  testing  was  completed,  field  measurements  were 
started  in  a cylinder  of  distilled  water  with  kh  = 1.5ir  and  ka  = 1.6, 
where  k is  the  wave  number  in  distilled  water,  h is  the  half-length  of 
the  cylinder,  and  a is  the  radius  of  the  cylinder.  These  measurements 
are  nearing  completion.  Measurements  in  salt-water  cylinders  are  then 
planned.  The  electrical  properties  (relative  permittivity  and  conductivity) 
of  the  salt-water  solutions  will  be  measured  using  a simple  coaxial-line 
setup  developed  previously  at  Harvard.^" 

The  theoretical  determination  of  the  electric  field  in  a finite  body 

when  illuminated  from  outside  is  difficult.  Some  attempts  at  numerical 

2 3 

solutions  have  been  made  in  recent  years.  ' These  and  other  possible 
methods  are  currently  being  investigated  and  will  be  applied  to  the  water 
cylinder  with  and  without  the  central  conductor.  It  is  hoped  that  the 
experimental  results  will  serve  as  the  basis  for  an  analytical -numerical 
solution  within  experimentally  determined  useful  ranges  of  the  parameters. 
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IV. 5 Experimental  Studies  of  Antennas  in  Plasmas.  M.J.  Miller,  Contract 

N00014-75-C-0648;  Research  Unit  11. 

The  performance  of  a cylindrical  antenna  of  half-length  h in  a 

material  medium  depends  on  both  8h  and  a/8.  In  a plasma  medium  a/8 

and  8 are  indirectly  related  to  the  ratio  of  the  operating  frequency  <d 

to  the  plasma  frequency  oj^  and  the  ratio  of  the  collision  frequency  V 

to  (Dp.  Changing  the  operating  frequency  of  a plasma-immersed  cylindrical 

antenna  within  a 10  to  1 range  can  change  8 more  than  3 orders  of 

magnitude  while  the  a/8  ratio  can  easily  vary  throughout  the  range 

0< a/8 <100.  A conference  paper  has  been  accepted  for  presentation^  which 

will  compare  the  results  of  an  experimental  investigation  to  measure  the 

impedances  and  distributions  of  current  and  charge  of  antennas  in  a plasma 

with  the  corresponding  theoretically  calculated  quantities.  Use  will  be 

made  of  the  general  theory  for  a cylindrical  antenna  in  a homogeneous, 

2 

isotropic,  dissipative  medium  which  applies  to  antennas  of  electrical 
length  0h  from  zero  to  more  than  2ir  immersed  in  a medium  whose  a/8 
ratio  can  vary  from  zero  to  10.  The  experimental  apparatus  used  to  measure 
the  electrical  performance  of  a plasma-immersed  antenna  consisted  of  a 
controlled  d.c.  helium  discharge  contained  in  a 15  cm  diameter  pyrex  tube. 
Plasma  frequencies  were  in  the  0.4-1  GHz  range,  and  v/to^  was  as  large 
as  0.2.  The  antenna  was  a bare,  11  cm  long,  cylindrical  monopole  driven 
from  a precision  coaxial  line.  It  was  fitted  with  small  movable  probes  to 
measure  distributions  of  current  and  charge.  The  coaxial  feed  line  had 
calibrated  voltage  and  current  probes  connected  to  a vector  voltmeter  which 
was  interfaced  with  a small  computer  so  that  antenna  impedances  could  be 
measured  and  plotted  during  experimentation. 
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IV. 6 Antennas  in  Matter:  Fundamentals,  Theory,  and  Applications.  R.W.P. 

King,  G.S.  Smith,  M.  Owens,  and  T.T.  Wu,  Contracts  N00014-75-C-0648 , 

F19628-76-C-0057 , and  F44620-72-C-0021,  and  NSF  Grant  ENG  75-14455; 

Research  Unit  11. 

The  results  of  current  and  past  researches,  supported  by  the  above 
Contracts  and  Grant,  on  bare  and  insulated  antennas  in  and  over  dissipative 
or  dielectric  media  have  been  integrated  into  a three-part  book  manuscript.1 
Part  I is  devoted  primarily  to  the  simple,  transmission-line-like  behavior 
of  the  insulated  antenna.  A brief  introduction  is  also  given  to  the  bare 
antenna  in  an  approximate,  quite  simple  form.  Both  types  of  antenna  are 
studied  as  probes  in  material  media.  The  propagation  of  plane  waves  and 
pulses  in  such  media  is  also  treated.  Part  II  is  an  advanced  theoretical 
treatment  of  antennas  in  various  media.  It  begins  with  a detailed  study 
of  the  electromagnetic  equations  in  homogeneous,  isotropic  media  including 
the  careful  formulation  of  constitutive  relations.  The  complete  theory  of 
the  bare  dipole  in  a general  medium  is  presented  and  this  is  followed  by  a 
rigorous  analysis  of  the  insulated  antenna.  Similar  analyses  are  provided 
for  the  bare  and  insulated  loop  antennas.  The  concluding  chapter  in 
Part  II  treats  the  problem  of  subsurface  communication  between  antennas 
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near  the  plane  interface  between  air  and  a general  half-space.  Finally, 
Part  III  is  concerned  with  the  construction  of  experimental  models  and  the 
techniques  of  measurement  relating  to  antennas  and  probes  in  a general 
dissipative  or  dielectric  medium. 
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IV. 7 Single  and  Crossed  Electrically  Short  and  Thin  Tubular  Cylinders. 

R.W.P.  King,  D.J.  Blejer,  and  B.H.  Sandler,  Contract  N00014-75-C-0648 
and  Grant  AFOSR-76-3073;  Research  Unit  11. 

The  surface  densities  of  current  and  charge  induced  on  tubular, 

highly  conducting  cylinders  and  intersecting  cylinders  by  a normally 

incident  plane  wave  depend  in  both  their  axial  and  transverse  distributions 

on  the  length  and  radius  of  the  individual  conductors.  When  ka<0.1,  the 

general  series  of  transverse  Fourier  components  is  well  approximated  by  two 

terms  for  both  E-  and  H-polarizations . When  kh>l,  ka<0.1,  the  zero- 

order  axial  current  dominates  and  thin-cylinder  theory  and  junction 

conditions  are  adequate  to  determine  the  axial  distributions.  A paper  has 
1 2 

been  written  and  presented  in  which  it  is  shown  that,  as  kh  is  reduced 
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to  satisfy  the  conditions  ka  < kh  « 1 , the  magnitude  of  the  rotationally 
symmetric  axial  current  decreases  to  small  values  and  the  first-order  axial 
currents  together  with  the  transverse  current  dominate.  All  of  these  have 
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nonrotationally  symmetric  transverse  distributions.  They  are  determined 
by  the  local  incident  magnetic  field  and  are  in  a large  measure  independent 
of  junction  conditions.  The  paper  points  out  the  significance  of  the 
surface  currents  and  charges  on  aircraft  illuminated  by  an  incident 
electromagnetic  pulse  containing  low-frequency  components  for  which  the 
dimensions  of  the  aircraft  cure  electrically  short.  For  them,  calculations 
of  currents  and  charges  by  thin-cylinder  approximations  and  junction 
conditions  are  inadequate  and  must  be  supplemented  with  or  even  replaced 
by  currents  and  charges  calculated  from  the  first-order  term  in  the 
transverse  Fourier  representation. 


References 

1.  R.W.P.  King,  D.J.  Blejer,  and  B.H.  Sandler,  "Current  Induced  on  Single 
and  Crossed  Electrically  Short  and  Thin  Tubular  Cylinders  by  a 
Normally  Incident,  Plane  Electromagnetic  Wave,"  to  appear  in  IEEE  Trans, 
on  Antennas  and  Propagation. 

2.  R.W.P.  King,  D.J.  Blejer,  and  B.H.  Sandler,  "Current  Induced  on  Single 
and  Crossed  Electrically  Short  and  Thin  Tubes  by  a Normally  Incident, 
Plane  Electromagnetic  Wave,"  presented  at  the  1978  USNC/URSI  Meeting, 
held  at  the  University  of  Maryland,  College  Park,  May  15-19,  1978. 


IV. 8 Surface  Currents  and  Charges  Induced  on  Electrically  Thick  Cylinders, 
Crossed  Cylinders,  and  Cylinders  Crossed  with  Horizontal  Flat  Plates 
by  an  Incident  Plane-Wave  Field.  R.W.P.  King,  D.J.  Blejer,  and  M. 

^ Krook,  Contract  N00014-75-C-0648  and  Grant  AFOSR-76-3073;  Research 

Unit  11. 

m 

A systematic  study  of  the  distributions  of  current  and  charge  induced 

on  crossed  metal  structures  was  begun  with  electrically  thin  wires  for  which 

12  3 4 

complete  analytical  ' and  experimental  ' solutions  were  obtained.  The 
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study  was  continued  with  crossed  electrically  thick  tubular  cylinders  for 

which  no  theoretical  analysis  is  available  so  that  a primarily  experimental 

investigation  was  carried  out. 6 However,  theoretical  distributions  of  the 

currents  and  charges  on  isolated  tubular  cylinders  with  arbitrary  lengths 

and  radii  were  evaluated6'7  with  E-  and  H-polarized  illuminations  both  for 

0 

standardizing  the  techniques  of  measurement  and  for  interpreting  the 
experimental  observations  on  crossed  cylinders.  The  initial  measurements 
were  carried  out  using  cylinders  and  crossed  cylinders  with  ka  = 1 . More 
recently,  these  have  been  extended  to  cylinders  and  crossed  cylinders  with 
ka  = 2 . 9 

The  next  major  step  in  an  approach  to  a model  of  an  aircraft  is  a 
study  of  the  currents  and  charges  on  the  surfaces  of  am  electrically  thick 
circular  cylinder  with  a transverse  member  consisting  of  a flat  plate.10 
With  the  aid  of  graphs  a detailed  description  is  given  of  the  experimentally 
determined  distributions  of  the  currents  and  charges  induced  by  a normally 
incident  plane  electromagnetic  wave  on  all  of  the  surfaces  of  a metal 
cross  formed  by  a vertical  tubular  cylinder  and  a horizontal  flat  plate. 

The  cylinder  is  a wavelength  in  circumference  and  1.75A  high  over  a large 
ground  plane.  The  horizontal  plate  consists  of  two  equal  arms  that  extend 
0.75A  from  the  surface  of  the  cylinder;  it  is  0.2125A  wide  and  0.0255X 
thick.  It  is  centered  at  a height  of  1.25A  from  the  ground  plane.  The 
distributions  are  compared  with  the  corresponding  ones  on  crossed  tub  iar 
cylinders  and  on  isolated  cylinders.  In  general,  the  currents  and  charges 
on  a horizontal  flat  plate  are  less  closely  coupled  to  the  currents  and 
charges  on  the  vertical  cylinder  than  are  those  on  a horizontal  cylinder. 

The  distributions  are  relatively  simple,  and  are  a useful  approximation  of 
the  surface  currents  and  charges  induced  on  an  aircraft. 
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IV. 9 Currents  and  Charges  on  Cylinders  in  a Parallel -Plate  Transmission 

Line.  R.W.P.  King,  D.J.  Blejer,  J.E.  Baum,  T.K.  Sarkar,  and  M.  Krook, 

Grant  AFOSR-76-3073  and  Contract  N00014-75-C-0648;  Research  Unit  11. 

In  order  to  determine  the  magnitudes  and  distributions  of  currents 
and  charges  induced  on  the  surfaces  of  an  aircraft  when  it  encounters  an 
electromagnetic  pulse  generated  by  a nuclear  explosion,  large  simulators 
have  been  constructed  in  which  the  aircraft  can  be  exposed  to  a hopefully 
comparable  electromagnetic  environment.  One  important  type  consists  of  a 
parallel-plate  transmission  line  or  open  waveguide  with  tapered  triangular 
ends,  between  the  conductors  of  which  a pulse  with  a plane  wave  front  must 
travel.  Since  the  relevant  pulse  contains  a frequency  spectrum  with  wave- 
lengths that  extend  from  very  long  to  relatively  short  compared  to  the 
dimensions  of  the  aircraft,  the  separation  of  the  two  conductors  of  the 
transmission-line  simulator  cannot  be  kept  small  compared  to  all  of  the 
wavelengths  in  the  pulse  as  required  for  TEM  transmission  without 
propagating  higher  modes.  It  follows  that  TEM-mode  transmission-line 
theory  can  be  applied  only  to  the  very  low  frequency  part  of  the  spectrum. 
The  properties  of  the  simulator  at  higher  frequencies  are  difficult  to 
determine  accurately.  One  approach  is  to  measure  the  electric  field  in  a 
model  simulator  at  one  or  more  of  the  relevant  higher  frequencies,  and  then 
seek  to  interpret  its  distribution  theoretically.  This  has  been  done  in  a 
carefully  designed  model  simulator  at  a frequency  of  625  MHz.*  3 When,  as 
in  this  case,  the  frequency  is  so  high  that  kh  is  not  small  compared  to 
one,  most  of  the  energy  supplied  by  the  generator  is  radiated  into  the 
surrounding  space,  and  only  a small  fraction  is  dissipated  in  the  load. 

The  structure  is  essentially  an  antenna  and  not  a transmission  line. 


Direct  measurements  were  made  of  both  the  amplitude  and  phase  of  the 

vertical  (E  ) and  longitudinal  (E  ) components  of  the  electric  field 
z y 

throughout  the  working  volume  (region  between  the  parallel  plates) , and  of 
Ez  on  the  surfaces  of  both  conducting  plates  where  E^=0.  In  order  to 
interpret  the  measured  results,  the  field  in  the  working  volume  was 
represented  approximately  in  terms  of  the  mode  theory  of  the  infinite 
parallel -plate  waveguide.  For  the  operating  wavelength  X = 48  cm,  four 
modes  cure  well  above  cut-off.  They  are  the  TEM,  ™02  an<^  ™03 

modes.  Comparisons  with  the  measured  data  indicate  that  the  assumed 
approximate  representation  of  the  field  by  the  superposition  of  ideal 
parallel -pi ate  modes  is  a good  one  throughout  the  working  volume  of  the 
simulator.  It  is  also  evident  that  the  total  E^  differs  greatly  from 
the  part  contributed  by  the  TEM  mode  alone.  This  is  true  particularly  near 
the  top  where  the  difference  between  the  incident  spherical  wave  in  the 
sloping  sections  and  the  plane  wave  in  the  working  volume  is  greatest. 

There  is  no  difference  on  the  ground  plane.  Note  also  that  the  field  in 
the  working  volume  is  not  a pure  traveling  wave  with  constant  amplitude, 
but  inclpdgs  a large  standing  wave.  TE  modes  contribute  only  negligibly 
jxid  are  ignored. 

Measurements  are  currently  being  taken  at  a second  high  frequency, 

f = 264  MHz,  where  with  the  separation  h between  the  parallel  plates 

reduced  from  108  cm  to  75  cm,  h is  now  .65X  rather  than  2.25X  of  the 

previous  case.  Now  only  the  TEM  and  TMQ1  modes  are  below  cut-off.  These 

can  be  uniquely  determined  from  the  longitudinal  and  transverse  components 

of  the  electric  field.  TE  modes  are  neglected.  The  TEM  mode  is  determined 

at  z/h=0.5  where  the  vertical  component,  E , of  the  electric  field  of 

z 


the  TH  , mode  is  zero.  E for  the  TM„,  mode  is  the  difference  between  the 

ox  z 01 

total  measured  E and  E for  the  TEM  mode.  The  measured  longitudinal 
z z 

component  of  the  electric  field,  Ey,  belongs  entirely  to  the  TMQ1  mode. 

An  attempt  was  made  to  measure  the  axial  component,  H^,  of  the  magnetic 
field  in  order  to  determine  the  magnitude  of  possible  TE  modes.  This  is 
difficult  since  a loop  antenna  must  be  used  as  the  probe.  Any  asymmetries 
in  the  loop  or  nonuniformities  in  the  field  due  to  the  presence  of  the 
edges  of  the  upper  plate  will  allow  unbalanced  currents  to  be  excited. 

Several  shielded  and  unshielded  loop  probes  have  been  constructed  and 
tested.  The  symmetry  of  the  probes  was  tested  by  comparing  the  response  of 
the  probe  in  a fixed  position  when  reversed  and  unreversed.  The  asymmetry 
of  the  probes  could  not  be  reduced  below  10%.  The  loop  probes  were  also 
tested  near  the  ground  plane  and  the  parallel  plate  for  their  ability  to 
detect  the  null  in  the  vertical  component  of  the  magnetic  field  required 
by  the  boundary  conditions.  A deep  null  was  detected  at  the  ground  plane 
but  not  at  the  parallel  plate.  It  appears  that  currents  induced  on  the 
conducting  tube  that  supports  the  loop  combined  with  small  asymmetries  in 
the  loop  make  it  impossible  to  measure  the  very  small  longitudinal  component, 
Hy,  of  the  magnetic  field. 

Another  area  of  investigation  during  this  reporting  period  has  been 
the  suppression  of  undesired  components  of  the  field  in  the  working  volume. 

As  a preliminary  representation  of  the  suggested  spatial  modal  filter,  a 
parallel  array  of  center-loaded  half-wave  dipoles  was  positioned  longitu- 
dinally in  the  working  volume  in  the  z = h/2  plane  for  different  values  of 
y.  No  attenuation  of  the  component  of  the  field  was  detected.  Further 

studies  are  to  be  carried  out. 
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IV. 10  Theoretical  and  Experimental  Study  of  the  Ferrite  Rod  Antenna. 


D.V.  Giri  and  R.W.P.  King,  Contracts  H00014-67-A-0298-0005  and 
N00014-75-C-0648,  and  the  Microwave  Physics  Branch,  Ballistic  Research 
Laboratories,  U.S.  Army  Aberdeen  Proving  Ground,  Maryland;  Research 
Unit  11. 

This  research  project  was  completed  some  time  ago.  During  this 
reporting  period  a technical  report  was  issued^-  which  presents  the  con- 
cluding part  of  the  study,  namely,  the  theoretical  analysis  of  practical 

2 

ferrite  rods  of  finite  length.  The  earlier  work  had  developed  a theory 
to  treat  the  electrically  small  loop  antenna  loaded  by  an  infinitely  long, 
homogeneous,  isotropic,  but  lossy,  ferrite  rod.  In  the  second  report,  two 
different  approaches  are  used.  The  first  method  makes  use  of  the  analogy 
between  the  ferrite  rod  antenna  and  the  conducting  cylindrical  dipole 
antenna  to  derive  an  integral  equation  for  the  magnetic  current  on  an 
infinitely  permeable  ferrite  antenna.  By  letting  h-*00,  the  integral 
equation  is  shown  to  be  consistent  with  the  earlier  results  for  the  infinite 
antenna.  Appropriate  corrections  to  the  integral  equation  are  then  made  by 


treating  the  ferrite  as  an  imperfect  magnetic  conductor,  and  an  approximate, 
three-term  expression  for  the  current  is  obtained.  The  second,  more 
rigorous  method  obtains  two  coupled  integral  equations  in  terms  of  the 
tangential  electric  field  and  the  tangential  electric  surface  current  from 
independent  treatments  of  the  interior  (ferrite)  and  exterior  (free  space) 
problems.  The  coupled  equations  are  then  solved  numerically  by  the  moment 
method.  Finally,  the  results  of  the  two  theories  are  compared  with  experi- 
mental measurements  made  on  eleven  different  antenna  configurations.  The 
agreement  is  good  for  all  cases  with  discrepancies  occurring  only  very 
close  to  the  driving  point. 
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